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Abstract. Properties of non-rotating and rapidly rotat- 
ing protoneutron stars and neutron stars are investigated. 
Protoneutron stars are hot, lepton rich neutron stars 
which are formed in Type-II supernovae. The hot dense 
matter is described by a realistic equation of state which 
is obtained by extending a recent approach of Myers and 
Swi§,tecki to the nuclear mass formula. We investigate the 
properties of protoneutron stars and neutron stars at dif- 
ferent evolutionary stages in order to emphasize the dif- 
ferences between very young and old neutron stars. The 
numerical calculations are performed by means of an exact 
description of rapid, uniform rotation in the framework of 
general relativity. We show that the minimal marginally 
stable protoneutron star mass is much higher than the 
corresponding minimum mass of a cold neutron star. The 
minimum gravitational (baryonic) mass of 0.89 - 1.13 Mq 
(0.95 - 1.29 Mq) of a neutron star is therefore determined 
at the earliest stages of its evolution. We also show that the 
use of different temperature profiles in the envelope as well 
as different shapes of the neutrino sphere change the prop- 
erties of protoneutron stars and hot neutron stars by up to 



1. Introduction 

It is generally believed that a neutron star (NS) is born 
as a result of the gravitational collapse of the iron core of 
a massive evolved progenitor star {M ^8 — 25 Mq) in a 
Type-II supernova (e.g. Bethe 1990). The iron core of such 
a star collapses when its mass reaches the Chandrasekhar 
limit 



Ch 



5.8Y^Mq - 1.2...1.5M, 



■ 



(1) 



where denotes the number of electrons per baryon 
which depends on the mass of the progenitor star. Due to 
the Fermi pressure of the nucleons, the collapse stops when 
nuclear matter density is reached and the core bounces 
back. Shortly after core bounce (some 10 ms) a hot, lep- 
ton rich NS, called protoneutron star (PNS), is formed. 
This PNS consists of a shocked envelope with an entropy 
per baryonpl s ~ 4 - 10 and an unshocked core with s ~ 1 



20%. I preliminary analysis indicates that even the most 
massive protoneutron stars rotating with Kepler frequency 
are secularly stable. Under the assumption of conserved 
angular momentum and baryonic mass, the maximum ro- 
tational frequency of an evolved neutron star is deter- 
mined by the Kepler frequency of the protoneutron star. 
We ca n thus derive a lower limit. P„,i„ ^ 1.56 — 2.22 ms. 



(Burrows et al. 1995). The envelope and the core contain 
nearly the same mass of about 0.6 - 0.8Mq (slightly de- 
pending on the mass of the progenitor star, see Burrows 
et al. 1995; Keil et al. 1996| ). During the so-called Kelvin- 
Helmholtz cooling phase (e.g. [lanka 19*93 ) the lepton num- 



ber decreases in the PNS due to the loss of neutrinos and 
consequently the PNS evolves in about 10 to 30 seconds 
into a hot, lepton poor neutron star (HNS) with an en- 
tropy per baryon, s ^ 1 — 2, depending on the model (e.g. 
Burrows & Lattimer 198(: ; Keil fc Janka 1995| ; ^umiyosh^ 



to the I rotational period of voung neutron stars with a 
canonical gravitational mass of 1.35 Mq. This result fur- 
therly supports the assumption that millisecond pulsars 
are accelerated due to accretion onto a cold neutron star. 

Key words ; Stars ; evolution - Stars ; neutron - Stars 



et al. 1995; Pons et al. 1999). After several minutes this 



HNS cools to a cold neutron star (CNS) with temperatures 



T < 1 MeV throughout the star (e.g. Keil fc Janka 1995) 



pulsars - Stars; rotation - Dense matter - Equation of — al. I99t ) 



Finally, it slowly cools via neutrino and photon emission 
until its thermal radiation is too weak to be observable 
after about 10^ yr (e.g. [Tsuruta fc Cameron 1966 ; Bchaat 



state 



Send offprint requests to: K. Strobel, 

e-mail; kstrobel@laser.physik.uni-muenchen.de 



The PNS is in /3-equilibrium during its lifetime, since 
the time scale of the weak-interaction is much smaller 
than the evolutionary time scale, i.e. the neutrino diffusion 
time scale or the neutrino cooling time scale, respectively. 
Hence the evolution of a PNS can be studied by consider- 
ing quasi-stationary models at different times. 



1 throughout this paper. 
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The properties of PNS's were investigated by differ- 
ent authors. For example, the case of no n -rotating PNS's 
was studied by Takatsuka et al. (1994 ), Bombaci et al 
(1995|3 and [Prakash et al. (19971). The cas e of rotating 



and the positron number density, respectively) and the net 
electron neutrino fraction, Y^^ — (ri^^, — np^)/n, (where n^^ 
and rip^ are the electron neutrino number density and the 
electron anti- neutrino number density, respectively). Since 
the muon number density is small in a PNS, they are ne- 
glected here. The reason for the small muon number is that 
the Kepler frequen c y which was developed for C NS's (e.g. the muon lepton family number is conserved, Y/^^+Y^^ = 0, 

while the neutrinos are trapped^ Relativistic calculations 
lead to the conclusion that NS's and PNS's are composed 
not only of nucleons and leptons but also of hyperons and, 



PNS's was treate d by Romero et al. (1992 ) and Takat- 
suka ( f995| |1996D by means of an empirical formula for 



Haensel et al. 1995). Hashimoto et al. (1994) and Goussard 
et al. jl997 ) account for rapid rotation by using an ex- 



act, general relativistic approach. Finally, Goussard et al 
(1998|)]have performed the case of differential rotation of 
PNS's. Most of these authors did not utilize an equation 



possibly, of nucleon isobars (see, e.g. Pandharipande 1971 
pchaffner fc Mishustin 1996| ; palberg fc Gal 1997|; |Hubei 



of stat |3 (EOS) of hot matter throughout the whole star et al. 1998 ). Nevertheless, we shall not take these particle 
(except Romero et al. 1992; Goussard et al. 1997, 1998| ), species into account. In view of the rather large uncertain- 
but used an EOS of cold matter for the envelope of the 
star instead. As we will show, this simplification leads to 
radii (Kepler frequencies) which are too small (large). 

The aim of this work is to study the properties of non- 
rotating and rapidly, uniformly rotating PNS's. We use 
an exact, general relativistic approach to rapid rotation 
(3chaab 1998). The hot dense matter is described by a re- 
cently devolped EOS ( Strobel et al. 199S| ), which is based 
on a modern parametrisation of the Thomas-Fermi ap- 
proach for finite nuclei and cold nuclear matter performed 



ties of the hyperon couplings, we shall, as a first approach 
to this problem, neglect these additional degrees of free- 
dom. In the following we select four different stages in the 
evolution of PNS's and NS's, namely at times ti ~ 50 — 100 
ms, t2 ~ 0.5 — 1 s, ^3 ~ 10 — 30 s, and t4 = some minutes 
after core bounce. 



by Myers fc Swiatecki (1990, 1991, 1996). We generalize 



this approach to hot dense matter by taking the thermal 
effects on both the kinetic and the interaction energy into 
account. In this way, we construct a set of EOS's with 
different profiles of the entropy per baryon and different 
lepton numbers. We can follow the evolution of the PNS 
into a CNS, by means of this set, at different evolutionary 
stages. We also investigate the influence of the location 
and the shape of the neutrino sphere as well as the influ- 
ence of the value of the temperature in the star's envelope. 

The paper is organized as follows. Firstly, in Sect. || 
we will briefly review the physics in the interior of PNS's 
and describe the different evolutionary stages of PNS's 
and NS's. Furthermore, we discuss the location and the 
shape of the neutrino sphere. The EOS's of PNS and NS 
matter are described in Sect. H, where we emphasize the 
influence of finite temperature and trapped lepton num- 
ber. The properties of rotating and non-rotating PNS's 
and NS's are presented in Sect. ^. Finally, discussion of 
our results and conclusions are given in Sect. |. 

2. Inside a protoneutron star 

A PNS differs in several respects from a CNS: At the be- 
ginning of its lifetime the PNS contains a high lepton num- 
ber, Yi, since the core is opaque with respect to neutrinos. 
A further difference is the high temperature which can- 
not be neglected with respect to the Fermi temperature 
throughout the whole PNS. We define, as usual, the lep- 
ton number, Yi = Y^. + Y^^ , as the sum of the net electron 
fraction, Yc ~ {n^- — nc+)/n, (where n, n^- , and no+ are 
the baryon number density, the electron number density. 



2.1. Protoneutron stars about 50 - 100 ms after core 
bounce 

This early type protoneutron star (EPNS) is characterized 
by a hot shocked envelope with an entropy per baryon of 
s ^ 4 - 5 for densities n < 0.02 fm~^, an unshocked core 
with s ~ 1 for densities n > 0.1 fm~^, and a transition 
region between these densities ( [Burrows et al. 1995 ). The 
entropy per baryon in the very outer layers of an EPNS 
is larger than s ~ 10. However, these layers have only a 
small influence on the EPNS structure and are therefore 
neglected. We investigate EPNS models with constant lep- 
ton number, Y\ = 0.4 for densities above n = 6x lO"** fm~^ 



where the neutrinos are trapped (Burrows et al. 1995). 
Below this density, the neutrinos can freely escape and 
the chemical potential of the neutrinos vanishes, fi^^ = 
(Coopcrstcin 1988). We refer to Table [| for the detailed 
parameters of the EPNS models studied here. 



2.8. Protoneutron stars at t 0.5 — 1 s after core bounce 

At this later stage, the entropy per baryon is approx- 
imately constant throughout the star, s ~ 2, except 



in some outer 


regions ( 


Burrows fc Lattimer 1986; 


Kei] 


|fc Janka 1995 




Keil 1996 


). The lepton number is ap- 



proximately constant since the neutrino diffusion time, 
Tiz-diff ~ 10 s, is by an order of magnitude larger than 
the PNS's age. We model this late type protoneutron star 
(LPNS) with a neutrino transparent envelope with densi- 
ties n < Ucnv = 6 X 10~^ fm~^ and a neutrino opaque core 



^ This statement is also true for massive neutrinos recently 
detected by the Super-Kamiokande collaboration (see Fukuds 



pt al. 1998| ). The detected mass is too small to allow consider- 
able flavour oscillations during th e first second s of a supernova 
(see discussion in Section 9.5.3 of Raffelt 199t) 
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Table 1. Entropies, temperatures, densities, and lepton numbers used in this paper. The entries are: entropy per baryon 
or temperature in the envelope, Scnv, Tonv; entropy per baryon or temperature in the core, Score, ?corc; maximum baryon 
number density of the envelope correlated with the entropy per baryon or temperature in the envelope, n(sonvj ?onv); 
minimum baryon number density of the core correlated with the entropy per baryon or temperature in the core, 
'^'(scorc, Tcorc); baryou number density below which the neutrinos are not trapped, n(ll,cnv); baryon number density 
above which the neutrinos are totally trapped, n(ll.coro); lepton fraction inside the core, 11, core- 



Label Scnv , TcYiv Score, TcoTC ^(^cnvi^cnv) ^(^corc , ^corc) ^(ll,cnv) ^(ll,corc) ll,corc 



[fn 



ffm" 



EPNS^£;,4 




5.0 


1.0 


0.02 


0.1 


0.0006 


0.0006 


0.4 


EPNS^to4 




4.0 


1.0 


0.02 


0.1 


0.0006 


0.0006 


0.4 


LPNSfL04 




2.0 


2.0 






0.0006 


0.0006 


0.4 






0.6 MeV 


2.0 


0.0004 


0.0004 


0.0006 


0.0006 


0.4 


LPNSyl04{64- 
LPNSyl04{64- 
LPNSyl04{64- 

LPNS^Los 


-63) 
-22) 
-62) 


2.0 
2.0 
2.0 
2.0 


2.0 
2.0 
2.0 
2.0 






0.0006 
0.0006 
0.0006 
0.0006 


0.006 
0.02 
0.06 

0.0006 


0.4 
0.4 
0.4 
0.3 


LPNS^Yo4 




1.0 


1.0 






0.0006 


0.0006 


0.4 


LPNS^Yos 




1.0 


1.0 






0.0006 


0.0006 


0.3 


LPNS$'lo4 














0.0006 


0.0006 


0.4 






2.0 


2.0 












jjNgT06s2 




0.6 MeV 


2.0 


0.0004 


0.0004 












MeV 


2.0 


0.06 


0.06 












1.0 


1.0 
















0.3 MeV 


1.0 


0.0004 


0.0004 












MeV 


1.0 


0.06 


0.06 








CNS 




MeV 


MeV 













* This notation means that the first number in parentheses classifies the lower density boundary of the neutrino sphere, i.e. 
6 X 10"*, the second number gives the upper density boundary, i.e. 6 x 10~^. 



with densities n > ricore and Yi = 0.4 (see Table |l|). The 
transition region between rienv and ricorc is called neutrino 



tion and shape of the neutrino sphere is rather uncertain. 
Another reason for the different choices is the uncertainty 



'T'core to simulate the influence of the shape of the neutrino coolin 


gofPNS's (e.g. 


Burrows & Lattimer 1988; 


Keil et al. 


sphere 


on the structure of LPNS's. We show also a LPNS 1996; 


Mezzacappa et al. 1998 


)■ 



with Yi — 0.3 for densities n > 6 x 10 fm for the sake 
of comparison, since different evolution calculations show 



different lepton numbers (Burrows & Lattimer 1986; Keil 



& Janla 1995; Pons et al. 1999) 



These different choices are motivated by recent theo- 
retical calculations of the neutrino-nucleon cross-section 
that include modifications due to the nucleon-nucleon in- 
teraction and spin-spin correlations. The problem was 
treated, for instance, by Raffelt (1996), Reddy et al. 
(1998|)] [Burrows fc Sawyer (1998D , and [Pons et al. (1999D 
using static correlation functions and by Raffelt (1996| ), 
[Janka et al. "(T996| ), ^Iffelt fc Strobel (1997| ), and T. Stro- 
bel (work in preparation) by using dynamical correlation 
functions. Because of the high complexity of the problem, 
the behaviour of the cross section and thus of the loca- 



For the sake of comparison, we investigate also a model 
with an isothermal envelope with temperature]^ T = 0.6 
MeV for densities below n ~ 4 x 10"'' fm~^. This temper- 
ature value is motivated by the fact that the temperature 
in the central parts of the progenitor star approximately 
raises to this value before the onset of the core collapse 
(e.g. [Shapiro fc Teukolsky 1983| ; [Bethe 1990[ ). ft is cer- 
tainly an upper limit of the true temperature since it cor- 
responds to an increase of the entropy per baryon by three 
or four orders of magnitude for densities n 6 x 10"^*^ 
fm^'^ (outer most layer of the star). This high entropy 



^ We choose T and not T* for simplicity, because the values 
of T and T* are comparable in the envelope; for the definition 



of T*, which includes metric corrections, see Soussard et al 



(1997) and Gondek et al. (1997) 
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seems to b e possible only in hot bubbles (see Mathews 
et al. l |993|) . 

For comparison we investigate also an unphysical, cold 
EOS, LPNSylo4: with a trapped lepton number of Y\ = 
0.4 for densities n > 6 x 10^* fm^^. 

2.3. Deleptonized hot neutron star at t ^ 10 — 30 s after 
core bounce 

After 10 — 30 s the neutrinos can freely escape and the 
HNS is nearly deleptonized. This also means that the 
lepton family number is not conserved anymore. The /3- 
equilibrium is thus given by /ip + /ic = Mn , Mc — f-fj.^ ^^id 
fi^ = for all neutrino species. At this stage muons have 
to be taken into account since the muon number density 
is comparable to the electron number density above nu- 
clear matter density. The entropy per baryon is nearly 
constant, s ^ 1 — 2, during the evolution from the LPNS 



to the HNS (e.g. Burrows & Lattimer 1986; 


Keil & Janka 


1995; 


Sumiyoshi et al. 1995; 


Pons et al. 19991). 



We again compare the models with isentropic en- 
velopes with models with isothermal envelopes, T = 0.3 
MeV or T 0.6 MeV for densities below n - 4 x 10""* 
fm~^. The models for the deleptonized HNS's are summa- 
rized in Table |l]. 



2.4- C^ld neutron star some minutes after core bounce 



After some minutes the NS has a temperature of T < 1 
MeV throughout the star and the EOS for cold NS mat- 
ter can be used to describe the CNS, because the ther- 



mal effects are negligibly small (see Shapiro fc Teukol 
sky 198^). W e shall adopt the model derived by Baym 



et al. (1971) for densities below neutron drip density, 
n< 2.6 X 10-4 and [Negele fc Vauthcrin (1973|) for 



densities between neutron drip density and the transition 
density, 2.6 x 10"'* < n < 0.1 fm^^^. Above this transition 
density, we use the model for CNS matter in /3-equ ilibrium 
without neutrinos derived by Strobel et al. (1997 ). 



3. Equation of state for protoneutron stars 

The EOS of PNS matter is the basic input quantity whose 
knowledge over a wide range of densities, ranging from 
the density of iron at the star's surface up to about eight 
times the density of normal nuclear matter reached in the 
cores of the most massive stars of a sequence, is neces- 
sary to solve the structure equations. Due to the high 
lepton number, the EOS of PNS's is different from the 
EOS's for cold and hot NS's with low lepton numbers. 
The nuclear EOS used in this paper for the description 
of a PNS is a Thomas-Fcrmi model of average nuclear 
properties, with a momentum- and density-dependent, ef- 



fective nucleon-nucleon interaction developed by Myers 



fe Swiatecki (1990, 1991). The parameters of the nuclear 
EOS were adjusted to reproduce a wide range of properties 



10 



10 



10' 



E 

> 10" 



10' 



10' 



EPNSv, 



CNS 



10" 



10"' 



10° 



10' 



10° 



10° 



Fig. 1. Pressure versus baryon number density for densi- 
ties n < 10~^ im~^ of hot dense matter. The curve CNS 



corresponds to cold matter (see Sect. 2A). The LPNSylo4 
and HNS""""®^^ curve is the isothermal envelope part of an 
isentropic core with entropy per baryon s — 2 (see Sects. 
U and U for explanation). Finally the curve EPNSf£^4 
corresponds to the envelope part of this EOS with s — 5 
and no trapped lepton number in this density region. 



of normal nuclear matter and nuclei (Myers & Swi§teck: 
1995| , |1996| , |1998| ; |Wang et al. 1997] ). ^trobel et al . (1999[ 



extended this approach to the case of finite temperature 
where they use exact numerical solutions for the inte- 
gration over the Fermi-Dirac distribution functions. Ap- 
pendix A contains a brief description of this approach. We 
extend the nuclear EOS to subnuclear densities and to dif- 
ferent compositions of PNS and HNS matter (i.e. trapped 
neutrinos, constant entropy per baryon, . . . ). In the sub- 
nuclear regime, the EOS was also obtained by means of the 
homogeneous Thomas-Fermi model. The electron number 
is derived by fitting the pressure to the subnuclear EOS's 
of |Baym et al. (1971] ) and [Negele fc Vautherin (197^ ) for 
densities below the density of the neutrino sphere in the 
case of EPNS and LPNS models and below the nuclear 
density in the case of HNS models. Our results for sub- 
nuclear densities are comparable to the EOS's derived by 



Lattimer fc Swesty (199l[) and used in the investigations 



of [Goussard et al. (1997| , |199^ ) and pondek et al. (1997] ). 

Figures |l| - ^ show the pressure, P(n), as function of 
the baryon number density, n, for different physical sce- 
narios. The envelope region with densities, n < 10~^ fm~^ 
is depicted in Fig. |l| for the isothermal part with T = 0.6 
of the LPNS^fof EOS and the HNS^''^'^ ^OS, the isen- 
tropic EPNS^5£^4 EOS, and the EOS for a CNS. The 
neutrinos do not contribute to the pressure for densi- 
ties n < 6 x 10~4 fm~^, since they are not trapped in 
this region. In the EOS with T = 0.6 MeV, the pressure 



^ The EOS's used in this paper are available at: http://www. 
physik . uni-muenchen .de /sektion / suessmann/astro . 
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10' 



10" 



10" 



E 

> 10"' 



lO'' 



10 



10 = 



EPNS" 




by a factor of ^ 3 — 4 each at n — 0.01 fm^'^. The ther- 
mal effects are even larger for the high entropy model la- 
beled EPNSylo4- The impact of the shape of the neutrino 
sphere can be inferred by comparison of the curves labeled 

LPNSylo4 ^-Ild LPNSylo4(64-22) ■ 

With increasing density the temperature dependence 
of the nuclear EOS increases. At the highest densities pos- 
sible in PNS's {n 1 — 1.2 fm^'^) the pressure increase 
due to thermal effects and due to high lepton numbers 
become comparable (see Fig. ^). This feature is clearly 
expressed by the nearly identical pressure of the HNS'^^ 
and the LPNS^t^ EOS's. 



n [ fm"^ 



The EOS of Lattimer & Swesty (1991) shows a smaller 
temperature dependence at high densities (see Goussar^ 



et al. 1997) in comparison with our EOS. The reason for 



this is that the temperature dependence of the EOS of 



Fig. 2. Pressure versus baryon number density in the den- 
sity region 10~^ fm~'^ < n < 10~^ fm~^ for different 
EOS's of hot dense matter. The abbreviations are de- 
scribed in Table |l|. 



1800 




Lattimer & Swesty (1991) lies entirely in the kinetic part 
of the enegy density, since they choose to* = m in their 
approach^. At this point it should be mentioned, that the 
behaviour of to* at high densities is highly uncertain. At 
n ~ 0.5 fm~^ we obtain, for instance, a pressure increase 
due to thermal effects of 28 % (from the LPNS^ t to 



-'YT-n4 

the LPNSylo4 case), whereas Lattimer fc Swesty (1991 ) 
obtain an i ncrease of only 8% (see Fig. l.b in Goussarq 
et al. 1997 ). The impact on the structure of PNS's and 
HNS's are discussed in Sec tion ^, where we compa re our 
results with the results of |Goussard et al. (1997, 1998) 



Swesty (1991) 



and Gondek et al. (1997) who used the EOS of Lattimei 



The number density and the mean energy of the neu- 
trinos y = Ue, i'e, V/^L, i^ii, i't, Vr are given bj|^: 

-dp 



9 
27r2 



P 



H-exp(i(p-/iy)) 



0.7 0.9 

n [ fm ' ] 

Fig. 3. Pressure versus baryon number density for densi- 
ties n > 0.1 fm~'^ for different EOS's of our model of hot 
dense matter. The abbreviations are described in Table |l|. 
The pressure of the HNS'^^ EOS is nearly identical to the 
LPNSylo4 case in this density region and is not shown for 
that reason. 



is dominated by the contribution of the photons, pph = 
2.817 X lO-'^T^MeV-^fm-s = 3.65 x lO-^MeVfm-^^ for 
low densities. Nevertheless, the influence of this low den- 
sity region on gross properties of LPNS's and HNS's is only 
small (see Sect. ^). In contrast to the isothermal EOS, the 
isentropic EOS is almost identical with the cold EOS in 
this density region. 

Figure || shows the general trend that the leptons dom- 
inate the pressure in the density region around and above 
the neutrino sphere, 10^* fm~'^ < n < 10~^ im^^ for 
LPNS's. Both thermal effects and the trapped lepton num- 
ber contribute significantly to the pressure and increase it 



and 



9 



2n,, 



p. 



(2) 



(3) 



la 1 + exp{^{p - fly)) 

respectively. The trapped electron neutrinos and anti- 
neutrinos are in chemical equilibrium, fip^ — —fiv^- The 
chemical potentials of all other neutrino species vanish, 
/ia; = (with X = Vfj^, Vfj^, Vr, i^r) duc to the lepton family 
number conservation (the muon number density is small 
compared to the electron number density in PNS's and is 
therefore neglected, for simplicity). The factor g denotes 
the spin-degeneracy factor and is related to the spin, s, of 
the particles by g = 2|s| -|- 1. Since only positive helicity 
neutrinos and negative helicity anti-neutrinos existFl the 



^ See Sect. 2.3 of Lattimer & Swesty (1991), espacially Eqs. 
2.8, 2.13, and 2.18. This will change if m* is chosen to be 
smaller than m. 



^ T he recently detected mass of neutrinos (see Fukuda et al 
1998| ) is much smaller than the mean energy of the neutrinos 



A finite neutrino mass therefore can be neglected. 

No negative helicity neutrinos nor positive helicity anti- 
neutrinos were found in experiment until now; CP violation of 
the weak interaction. 
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Fig. 4. Chemical potential versus baryon number density 
for electrons and electron-neutrinos of the LPNSylo4 ^^'^ 
LPNSylo4 cases. The dotted lines correspond to the CNS, 
the HNS'^i and the HNS''^ eOS's. 
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Fig. 6. Mean neutrino energies versus baryon number den- 
sity for all neutrino types in the case LPNSylo4 for den- 
sities larger than 6 x 10~* fm~'^ {ly^ = v^ji.,v^,Vr and Dt). 



> 






and 



"p, = — exp( — — ) 

IT- 1 



jexp — —). 



T 



(6) 



(7) 



This leads to the simple, linear temperature dependence 
of the mean energy of the electron anti-neutrinos: 



ST, 



(8) 



the energy of an ultra-relativistic Boltzmann gas. The 
mean neutrino energies of all neutrino species are shown 
in Fig. ^for the case LPNSylo4- Figure |4| shows the chem- 
ical potential of the electron neutrinos and the electrons 
for different models of the EOS. 

In Fig. H we show the density dependence of the tem- 
perature for different EOS's. One can see the temperature 
drop in the EPNSylo4 EOS at the interface between the 
hot shocked envelope {n < 0.02 fm^*^) and the unshocked 
degene racy factor is equal 1 for the neutrinos. In the case core (n > 0.1 fm ~^)(see [Burrows fc Lattimer 1986|; |Bur 



Fig. 5. Temperature versus baryon number density for dif- 
ferent EOS's. 



of vanishing chemical potential Eqs. (2) and (3) lead to a 
temperature dependence of the number density: 



rows et al. 1995 ). The temperature increases wit h inc reas- 

and 



4.2 



1.19 X 10" 



' MeV"^ fm 



-3 



(4) 



and a linear temperature dependence of the mean neutrino 
energy: 



3.15 T, 



(5) 



with X = Vfj^, D^, Vr, Vt. Due to the high chemical potential 
of the electron neutrinos in the case of a high trapped lep- 
ton number {^ly^ ^ T; see Fig. ^ and Fig. H), the number 
density and the mean energy of the electron anti- neutrinos 
can be approximated by: 



ing entropy per baryon (see discussion in Sect, 
decreases with increasing lepton number, since the neu- 
tron ft^^U^na^dthfiproton fraction become more equal 
(see Prakash et al. 1997| ). The maximum temperature in 
the most massive PNS's and HNS's reaches values between 
80 and 120 MeV (see Fig. | and Table |). The tempera- 
ture in PNS's and HNS's with a typical baryonic mass of 
1.5Mq has values between 20 and 40 MeV (see Table H). 

The fractions of electrons and electron-neutrinos are 
shown in Fig. | for the LPNS^1lo4 ^nd LPNS^\o4 cases 
with constant lepton number. The lepton fractions of the 
unphysical LPNSylo4 EOS are nearly identical to that 
of the LPNSylo4 EOS and are not shown in Fig. ^ for 
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Fig. 7. Lcpton fractions versus baryon number density 
for different EOS's. The lepton fractions of the unphys- 
ical LPNSylo4 EOS are nearly identical to that of the 
LPNSylq4 EOS and is not shown for that reason. The 
lower dotted line corresponds to the CNS EOS, the mid- 
dle dotted hne to the HNS'^i EOS and the upper dotted 
hne to the HNS''^ EOS. The dotted lines show the sum 
of the electron and muon fraction. The stars correspond 
to the positron fraction of the HNS*'^ case. The positron 
fraction of all other EOS's lies below the resolution of this 
figure. 



that reason. The fractions are nearly constant (29 — 32% 
for electrons and 8 — 11% for electron neutrinos) in a 
wide range of densities. The electron fraction decreases for 
n > 0.7 fm~'^ since t he symmetry energy decreases for den- 
sities n > 0.4 fm~^ dStrobel et al. 1997[ ). Non-relativistic 
EOS's derived within variational approaches behave sim- 
ilar in this respect, whereas the symmetry energy derived 
in relativistic and non-relativistic Briickner-Bethe calcu- 



lations monotonically increases with density (see Strobel 
et al. ;997). Higher temperatures cause a small decrease 
(increase) of the electron (electron neutrino) fraction (see 
[Takatsuka et al. 1994 ). Without trapped neutrinos the 
sum of the electron and muon fraction increases, in first 
approximation, quadratically with increasing temperature 
(see Fig. | and [Keil fc Janka 1995| ). 

The effective masses, m*/mT, of neutrons and pro- 
tons for symmetric and pure neutron matter are shown 
in Fig. ^. The effective mass for symmetric nuclear matter 
at saturation density is found to be m*/mr — 0.867 for 
neutrons and protons in our model. This is in good agree- 
ment with the generally accepted, experimental value (e.g. 
Bauer et al. 1982). In the case of pure neutron matter, the 
effective neutron mass increases up to m*/TOn = 0.935, at 
nuclear matter density, whereas the effective proton mass 
in pure neutron matter decreases to m*/mp — 0.808. 
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Fig. 8. Effective mass of neutrons and protons versus 
baryon number density. The solid line shows the effective 
mass of neutrons and protons for symmetric nuclear mat- 
ter (the difference in the values for neutrons and protons 
are negligibly small). The long dashed line shows the ef- 
fective mass of neutrons in pure neutron matter. The dot- 
dashed line shows the effective mass of protons brought 
into pure neutron matter. 




Fig. 9. Adiabatic index versus baryon number density for 
different EOS's. 



The adiabatic index 7: 



dlnP 



7 = 



din? 



n dP 
P d^ 



(9) 



is shown for different EOS's in Fig. ^j. The adiabatic index 
decreases with increasing temperature and lepton numbers 
for densities around and above nuclear density. In con- 
trast, it decreases with increasing lepton number for den- 
sities n < 0.05 fm~'^ (see Gondek et al. 1997). The steep 
behaviour of the adiabatic index 7 around nuclear matter 
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Table 2. The speed of sound Vs in the density region around and above nuclear matter density for different EOS's. 



The abbreviations are described in Table ^ The maximum value is reached in the CNS EOS: 
fm^^ Vs ~ 0.964 (in units of c). 



^(CNS) = 1.246 



n [fm 



EOS 


0.02 


0.06 


0.1 


0.2 


0.3 


0.4 


0.6 


0.8 


1.0 


1.2 


LPNSfL04 

LPNS|.Yo4 
CNS 


0.1242 
0.1046 
0.1011 
0.0684 
0.0507 


0.1627 
0.1320 
0.1473 
0.1095 
0.0807 


0.2067 
0.1737 
0.1917 
0.1559 
0.1255 


0.3270 
0.2968 
0.3128 
0.2735 
0.2469 


0.4411 
0.4172 
0.4280 
0.3905 
0.3699 


0.5434 
0.5240 
0.5305 
0.4962 
0.4802 


0.7089 
0.7008 
0.6963 
0.6665 
0.6569 


0.8276 
0.8257 
0.8176 
0.7921 
0.7863 


0.9055 
0.9113 
0.9073 
0.8862 
0.8812 


0.9594 
0.9617 
0.9750 
0.9606 
0.9543 



density is the reason for the core b ounce of the coUaps- 
ing iron core of the progenitor star ( ^hapiro fc Teukolsky 
1983()7] 

The speed of sound Vs in units of the speed of light c: 




<'»> 

is tabulated in Table || for different EOS's. The speed of 
sound Vs increases with density up to nearly the speed of 
light, in the most massive stars of a sequence. Neverthe- 
less, Vg is always smaller than the speed of light and all 
EOS's used in this paper are causal. The speed of sound 
increases with temperature and lepton number at fixed 
density in contrast to the results of Goussard et al. (1998] ) . 
This is probably caused by the smaller temperature depen- 



dence of the EOS of Lattimer & Swesty (1991 ), which was 
discussed before in this section. 

4. Structure of rotating and non-rotating 
protoneutron stars 

The structure of rotating PNS's and NS's is governed by 
the Einstein equations in stationary, axisymmetric, and 
asymptotic flat space-time. Under these special condi- 
tions to the space-time symmetry the ten Einstein equa- 
tions reduce to four non-trivial eq uations which are ellip - 
tic in quasi-isotropic coordinates ( Bonazzola et al. 1993| ). 
The four non trivial Einstein equations together with the 
energy-momentum conservation are solved via a finite dif- 
ference scheme (Schaab 1998). We follow Bonazzola et al. 



(1993 ) in compactifying the outer space to a finite region 
by using the transformation r — > 1/r. The boundary con- 
dition of approximate flatness can then be exactly fulfilled. 
The neutron star model is uniquely determined by fixing 
one of the parameters: central density, gravitational mass, 
or baryon number, as well as one of the parameters: an- 
gular velocity, angular momentum, or stability parameter 
p = £'kin/|-£'grav|- The modcls of maximum mass and/or 
maximum rotation velocity can also be calculated. 




ejgcm' ] 

Fig. 10. The gravitational mass versus central energy- 
density of non-rotating NS's and PNS's. The lower long 
dashed line corresponds to the HNS''^ EOS and the upper 
long dashed fine to the HNS^^ EOS. The lower solid line 
corresponds to the LPNSylo4 EOS and the upper solid 
line to the LPNS|?lo4 EOS. 



Protoneutron star and neutron star sequences 

Figure |l^ shows the gravitational mass as function of the 
central energy-density for the non-rotating PNS and NS 
models. Only models whose gravitational masses increase 
with central energy-density are stable against axisymmet- 
ric perturbations. Whereas the maximum mass differs only 
by 5 % between the various star models, the mini- 
mum mass of the EPNS models, Af,„i„ 0.9 - 1.2Mq, 
is much larger than the minimum mass of the CNS mod- 
els, M^in - O.IMq. 

Figure |l^ shows the Kepler frequency, i.e. the fre- 
quency at which mass shedding sets in, as function of the 
gravitational mass. Also shown is the rotational frequency. 
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Fig. 11. Kepler frequency versus gravitational mass of 
PNS's and NS's for different EOS's. The upper long 
dashed line corresponds to the HNS*^^ EOS and the lower 
long dashed line to the HNS^^ EOS. The upper solid hne 
corresponds to the LPNSylq4 EOS and the lower solid line 
to the LPNS^\o4 EOS. 
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Fig. 12. The gravitational mass versus central energy- 
density of NS's and PNS's rotating at their Kepler fre- 
quency. The lower long dashed line corresponds to the 
HNS'^i EOS and the upper long dashed line to the HNS^^ 
EOS. The lower solid line corresponds to the LPNSylq4 
EOS and the upper sohd line to the LPNS^\o4 EOS. 



1 7 = 4033 s~\ of th e fastest pulsar known, PSR 1937+214 
( Backer et al. 198^ ). The gravitational mass of this pulsar 
is unfortunately unknown but is typically assumed to be 
in the range 1.0-2.0 Mq. Since the Kepler frequency of 
the CNS models is larger than the rotational frequency of 
PSR 1937+214 for this mass range, this model is consis- 
tent with the observation. 

Figure |l^ shows the sequences of stars rotating with 
Kepler frequency. As expected the masses are now larger 
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Fig. 13. The gravitational mass versus stellar radius (as 
measured by an observer located at infinity) of non- 
rotating NS's and PNS's. The abbreviations for the differ- 
ent EOS's are described in Table |l|. The lower long dashed 
line corresponds to the HNS^^ EOS and the upper long 
dashed hne to the HNS**^ EOS. The lower solid hne corre- 
sponds to the LPNSylq4 EOS and the upper solid hne to 
the LPNS^\o4 EOS. 
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Fig. 14. The gravitational mass versus equatorial radius 
(as measured by an observer located at infinity) of NS's 
and PNS's rotating at their Kepler frequency. The lower 
long dashed line corresponds to the HNS*^^ EOS and the 
upper long dashed line to the HNS'^^ EOS. The lower solid 
line corresponds to the LPNSylo4 EOS and the upper 
sohd line to the LPNSfL04 EOS. 



at fixed central energy density, but the maximum mass 
is reached for lower central energy densities (see also Ta- 
ble |). The mass increase due to rotation is smaller for 
PNS's than for CNS's, since the radii of PNS's are larger 
and thus the Kepler frequencies are smaller. 
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Fig. 15. Iso-energy-density surfaces for a EPNSylo4 
model with a baryonic mass of 1.5Mq rotating 
at its Kepler frequency. The surfaces correspond 
to energy-densities (beginning at the centre), e — 
0.99, 0.1, 0.01, 10-3, 10-^ 10-5, lO-'^Ec, where = 8.07 x 
g cm" 



10^"* g cm-'^ denotes the central energy-density. 



We also show the mass-radius relation for PNS's in 
In comparison with CNS's, where the 



Figs, y and |14 



radius only slightly changes in the relevant mass region 
around Mq ~ 1.4 A/q, one obtains for PNS's a much 
stronger increase of the radius with decreasing mass, 
which is caused by the stiffer EOS's for the PNS's. 

To demonstrate the differences of the internal structure 
of an EPNS and a CNS, we show in Figs. H and |l| the 
iso-energy-density surfaces for models with Mb — I.SA/q. 
It turns out that EPNS's contain mainly matter of den- 



sities below nuclear matter density (Fig. 15). Contrary to 
this, matter in a CNS is dominated by matter of densities 
around and above nuclear matter density (Fig. |l^) . 

Another important point concerns the stability of the 
star models against secular or dynamical instabilities. A 
certain configuration is dynamically stable against ax- 
isymmetric perturbations if the gravitational mass is min- 
imum with respect to variations at fixed baryon num- 
ber and angular momentum. Along a star sequence with 
fixed angular momentum, this is the case if the gravi- 
tational mass increases with the central energy density 
(see Shapiro fc Teukolsky 1983| , p. 151). Please note that 
the angular momentum is not fixed in the sequences of 
stars rotating with Kepler frequency. The maximum mass 
configuration differs therefore slightly from the marginally 
stable configuration. A fully general relativistic analysis of 
dynamical and secular stability against non-axisymmetric 

The equatorial radius _Rinf is the circumferential radius as 
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Fig. 16. Iso-energy-density surfaces for a CNS 
model with a baryonic mass of 1.5A/q rotat- 
ing at its Kepler frequency. The surfaces corre- 
spond to energy-densities (beginning at the centre), 
£ = 0.99, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.l£c, where 
Ec = 8.11 X lO^'' g cm^3 denotes the central energy- 
density. 



perturbations is extremely difficult and has been per- 
formed only by means of approximations and/or special 
assumptions in literature. It turned out however that the 
configurations can be classified by a stability parameter 
defined by the ratio of the kinetical energy and the abso- 
lute value of the gravitational energy, fi = £'kin/|-Egrav|- 
For values, (3 > 0.26, the models are probably dynam- 
ically unstable against the bar mode (to = 2), whereas 
they become already secularly unstable for (3 > 0.14 (see, 
e.g., Durisen 1975; Managan 1985| ; [mamura et al. 1985; 



Bonazzola et al. 199(]| ). Higher modes 



measured in infinity (see, e.g., Schaab 1998) 



TO > 5, of pertur- 
bations are suppressed by the finite viscosity of neutron 
star matter. It is however uncertain whether modes with 
TO = 3 or 4 are also suppressed. If this is not the case, 
secular in stability sets already in above a critical value, 
l3 > 0.08 ([Friedman et al. 1986| ). 

Figure ^ shows the value of the stability parameter 
for stars rotating with their Kepler frequency. Since the 
timescale for the growing of the secular instabilities is 
larger than the evolution timescale of PNS's and HNS's, 
only the dynamical instability is important for PNS's and 
HNS's. The critical value for the dynamical instability, 
P > 0.26, is reached by none of the models. Depending 
on the internal structure of the CNS's, the CNS's models 
rotating with Kepler frequency may become secular un- 
stable against modes with to = 3 or to = 4 for star masses 
M > O.SMq. 
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Fig. 17. Stability parameter /3 = £'kin/|£'grav| versus grav- 
itational mass of PNS's and NS's rotating at their Ke- 
pler frequency for different EOS's. The upper long dashed 
line corresponds to the HNS'^^ EOS and the lower long 
dashed line to the HNS''^ EOS. The upper sohd line cor- 
responds to the LPNSylo4 EOS and the lower solid line 
to the LPNS^\o4 EOS. 



2.2 
2.0 
1.8 

"1 ■'■6 
"S 1.4 
1.2 
1.0 
0.8 



0.8 



CNS 
ERNSv 



1.0 1.2 1.4 1.6 1.. 

M„ [ M„,„ 1 



2.0 2.2 2.4 2.6 



Fig. 18. Gravitational mass versus baryonic mass for the 
non-rotating EPNS^£^4 EOS and the non-rotating CNS 
EOS. The EPNS^Ssi^ ^Qg j-^g exactly on the EPNS^^m 
line and is not shown by that reason. The arrows show 
the evolution (from left to right) for the minimum mass 
EPNS^£i4, the minimum mass EPNS^5£i4, the Mb = 
1.5 Mq, and the maximum mass configurations. 



Evolution of a non-rotating protoneutron star 

The ev [3hition of a PlNb to a CJiNS can be toUowed by means 
of several "snapshots" taicen at ditterent tunes after core 
bounce (see Sect. ||). The evolution path is determined 
by fixi ng the baryonic mass and the angular momentum 
if accretion of matter and loss of angular momentum is 
neglected. First, we study the evolution of a non- rotating 
PNS with Mb = 1.5M0 (see Table | and Fig. |l|). 

The first snapshot corresponds to approximately 50 
- 100 ms when the envelope is characterized by a high 
entropy per baryon and high lepton number (see curves 
labeled EPNSf£i4 and EPNS^^m)- After 0.5 - 1 seconds 
the PNS reached our LPNS stage, which is characterized 
by an approximately constant entropy per baryon, s 2, 
throughout the star (model LPNSylo4)- Due to the higher 
entropy per baryon in the core of the PNS, the central den- 
sity and thus the gravitational binding energy decreases. 
On the other hand the entropy per baryon decreases in 
the envelope and therefore the radius decreases, too. For a 
lower value of the entropy per baryon and/or a lower value 
for the lepton number, the gravitational binding energy in- 
creases compared to the EPNS stage (models LPNSylosj 
LPNS^Yo4' LPNS^^Loa)- ^^^r about 10 - 30 seconds the 
neutrinos escape from the star and its EOS softens (HNS 
models). The gravitational binding energy thus increases 
by roughly 3 %. Finally, the CNS model is even more com- 
pressed. 

Due to a smaller increase of the pressure with increas- 
ing temperature. Burrows & Lattimer (1986) and Pons 



tral density during the evolution of the PNS to the CNS. 
In contrast, our results, and also those of Keil & Janka 
(1995D , show a more complex behaviour of the central 



density (see discussion in the paper of Pons ct al. 1999) 



et al. (1999) obtain a monotonous increase of the cen- 



This difference has mainly two reasons: Firstly, Pons et al 
(1999D use an entropy profile obtained in a supernova col 
lapse simulation of a 1.08 Mq NS, which leads to a higher 
central entropy per baryon, s ~ 1.5, at the EPNS stage 
and an increase to 3 2 at the LPNS stage. We use 
the calculations of Burrows et al. (1995 ) who simulated 
the supernova collapse of a NS with M ~ 1.5 Mq. They 
obtained a smaller starting central entropy per baryon, 
s ~ 1, and an increase to s ~ 2 at the LPNS stage. The 
second reason is the fact, that they use an approximation 
for the temperature influence on the pressure (derived by 
Prakash et al. 1997) of the nucleons: 



T, 



(11) 



Po 



5s 



3 d Inm* 
2 d InriT- 



37r2 



pot 



kin 



(12) 



where T-p,T denotes the Fermi temperature of the quasi 
particles (7f,t = Pft/^'ti^), Pth, the thermal pressure 
and Po — Pkin + Ppot, the pressure at zero temperature 
as sum of kinetic and potential pressure. This approx- 
imation holds under the assumption that in the dense 
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Table 3. Properties of non- rotating and with Kepler frequency rotating PNS's and NS's, for a fixed baryonic mass 
Mb — 1.5 Mq. The EOS's are summarized in Table |^. The entries are: gravitational mass, Mg', baryonic mass, Mb; 
circumferential radius (as measured in infinity), i?inf; central baryon number density, tIc; central temperature, Tc] 
Kepler frequency, fix; angular momentum, J (Mq km =5.966 x 10''® g cm^ s~^). 



n = 



EOS 




Mg 


Rint 




Tc 


Mg 


-Rinf 






J 






[Mq] 


[km] 




[MeV] 


[Mq] 


[km] 


[fm-^] 




[Mq km] 






1.425 


25.66 


0.440 


20.8 


1.430 


38.55 


0.427 


1879 


0.634 






1.425 


18.66 


0.456 


21.3 


1.433 


27.46 


0.436 


3056 


0.924 


LPNS^io4 




1.431 


17.38 


0.385 


39.5 


1.437 


25.62 


0.353 


3410 


1.234 


LPNSp? 




1.432 


17.47 


0.386 


39.6 


1.437 


25.77 


0.355 


3381 


1.219 


LPNSyl04{64- 


-63) 


1.430 


16.72 


0.386 


39.6 


1.438 


24.68 


0.348 


3610 


1.340 


LPNSyl04(64- 


-22) 


1.430 


16.34 


0.386 


39.6 


1.438 


24.13 


0.344 


3739 


1.408 


LPNSyl04{64- 


-62) 


1.430 


15.85 


0.388 


39.7 


1.439 


23.36 


0.341 


3930 


1.497 


LPNS^Los 




1.414 


16.01 


0.400 


42.9 


1.423 


23.37 


0.361 


3897 


1.362 


LPNS^Yo4 




1.411 


14.70 


0.459 


21.4 


1.423 


21.26 


0.417 


4482 


1.389 


LPNS^Yo3 




1.391 


13.59 


0.475 


22.8 


1.406 


19.55 


0.424 


5076 


1.552 


LPNS^Olo4 




1.410 


12.39 


0.519 





1.429 


17.73 


0.450 


5949 


1.800 






1.391 


15.22 


0.412 


49.7 


1.403 


22.13 


0.367 


4204 


1.431 






1.392 


15.29 


0.412 


49.7 


1.404 


22.24 


0.369 


4171 


1.413 






1.385 


13.53 


0.404 


49.0 


1.407 


19.86 


0.324 


5075 


2.072 






1.362 


12.78 


0.501 


26.8 


1.382 


18.36 


0.437 


5562 


1.669 






1.363 


12.80 


0.502 


26.9 


1.382 


18.39 


0.438 


5550 


1.662 






1.362 


12.18 


0.504 


27.0 


1.387 


17.59 


0.418 


6023 


1.957 


CNS 




1.346 


11.45 


0.551 





1.374 


16.36 


0.459 


6689 


1.953 



parts of the PNS the temperature is small in compari- 
son with the Fermi temperature, i.e. T/Tp <^ 1. However 
in our case these ratio reaches T/Tp w 0.45 {T/T-p k, 0.6) 
at n = no and T/Tp « 0.29 (T/Tp « 0.44) at n = 1 
fm^"^ for neutrons (protons) in the LPNSylo4 mo del (s ee 
also the Fermi-Dir ac distribution function in Fig. A.l of 
Appendix A and Takatsuka et al. (1994 ) for this pur- 
pose, in which it can be seen that more than 10 % of 
the matter is non-degenerate). Therefore this approxima- 
tion underestimates the pressure increase due to thermal 
effects, e.g. the thermal pressure Pth is exact (approxi- 
mative) 2.37 MeVfm-3 (1.26 MeVfm'^) at n = no and 
54.62 MeVfm-3 (28.83 MeVfm'^) at n = 1 fm'^ for the 
nucleons of the LPNSylo4 '^^^ of our model. Hence in 
our opinion, higher order terms should be included in the 
treatment. Since the temperature dependence of the GM3 
model of Pons et al. (1999) is smaller than in our case, the 
deviations between the exact solution and the approxima- 
tion may be small. 



^.3. M hximum rotational frequency of a neutron star 



state. During its evolution, the angular momentum is as- 
sumed to be conserved (see Table As in the non- 
rotating case, the star becomes more and more compressed 
during its evolution. With only a few exceptions, the cen- 
tral baryon density increases, whereas the gravitational 
mass and the circumferential radius decreases. It is ob- 
vious that this trend has to be counterbalanced by an 
increasing angular velocity in order to keep the angular 
momentum constant. Compared to the EPNS state, the 
angular velocity in the CNS state is increased by 51 % or 
32 % for the sequences starting with the model EPNSylo4 
or EPNSylo4i respectively. Nevertheless, the angular ve- 
locity in the CNS state reaches at most 60 % of the Kepler 
frequency. 

In this respect, the minimum rotational period is de- 
termined by the Kepler rotating EPNS model. Figure [l^ 
shows the the general evolution of the rotational period for 
the two EPNS models. For a CNS with a typical baryonic 
mass of 1.5Mq one obtains a minimum rotational period 
P, be t ween 1.56 and 2.22 ms. Recently, Goussard et al 



We follow now the evolution of an EPNS with Mb — 
1.5M0, which rotates with Kepler frequency at its EPNS 



(1997, 1998) found similar results for the minimum rota- 
tional period, which confirms our calculations. The period 
of the fastest known pulsar PSR 1937+214, P = 1.56 ms 



(Backer et al. 1982) is at the lower limit of this range. Pul- 
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Table 4. Evolution of 1.5 Mq star at constant angular momentum, J = 0.634 Mq km and J 
sequences starting with the model EPNSylo4 ^^^d EPNSylo4j respectively. 



0.924 Mq km for the 



EPNSlFfL 



EOS 



Mg 

[Mq] 



[km] 



71c 



[MoV] 



Mg 
[Mq] 



Rini 

Ikm] 



ric 
[fm" 



[McV] 



EPNS 


1.430 


38.55 


0.427 


20.4 


1879 


1.433 


27.46 


0.436 


20.7 


3056 


LPNSfL04 


1.349 


14.83 


0.350 


37.0 


2091 


1.341 


15.40 


0.336 


36.0 


2961 


LPNS^Los 


1.363 


14.27 


0.374 


40.9 


2171 


1.361 


14.74 


0.362 


39.9 


3073 


LPNS^Yo4 


1.391 


13.17 


0.442 


20.9 


2377 


1.391 


13.54 


0.430 


20.5 


3390 


LPNS^Yos 


1.385 


12.76 


0.462 


22.4 


2474 


1.387 


13.09 


0.451 


22.0 


3525 




1.349 


14.06 


0.387 


46.5 


2246 


1.348 


14.53 


0.375 


46.3 


3178 




1.362 


12.49 


0.489 


26.3 


2609 


1.365 


12.83 


0.478 


25.9 


3709 


CNS 


1.348 


11.83 


0.535 





2833 


1.351 


12.14 


0.523 





4031 



CNS(n,) 
EPNS=;,;;(a,) 

CNS(EPNS^t",(aj) 

EPNSt„;(a,) 

CNS(EPNS;l;,(a,)) 




Fig. 19. Minimum rotational period P versus baryonic 
mass of the EPNS's and CNS's. The figure shows the 
spin up during the evolution of the EPNS's to CNS's. 
The dot-dashed (short dashed) line shows the EPNSy£q4 
(EPNSylo4) model rotating at its Kepler frequency. The 
solid (long dashed) line shows the corresponding CNS 
model obtained by cooling at constant angular momen- 
tum. The dotted line shows CNS's rotating at their Ke- 
pler frequency for comparison. The arrows represent the 
evolution from the EPNS stage to the CNS stage. 



sars that rotate even faster cannot be born with such small 
periods but have to be accelerated after their formation, 
as long as a typical baryonic mass, Mb — 1.5Mq is as- 
sumed. Andersson (1998 ) has recently shown that not too 
cold NS's rotating with fl > 0.1 I^k are unstable against 
r-modes. This means that the minimum rotational period 
for a young NS is even higher, Pmin ^ 10 ms. Both the re- 



0.08 

0.07 
0.06 
0.05 
CO. 0.04 
0.03 
0.02 
0.01 



EPNS™„(n,: 
--- CNS (EPNs;:„(nK)) 

EPNSt>,: 

— CNS (EPNs;:„(n,)) 



0.9 



1.1 




Fig. 20. Variation of the stability parameter /3 = 
^^kin/l-Bgravl during the evolution from EPNS's to CNS's 
at constant angular momentum. The arrows represent the 
evolution from the EPNS stage to the CNS stage. 



strengthen the vie w of milliseco nd pulsars as being recy- 
cled by accretion ( Lorimer 1996 ). 

Figure ^ shows the behaviour of /3 during the evolu- 
tion from the EPNS to the CNS state. The stability pa- 
rameter increases during the evolution. Nevertheless, its 
maximum value, /3 ~ 0.082, reached for the most massive 
star of the EPNS^£^4 EOS, is only a little larger than 
the critical value f3 ^ 0.08 for the onset of the secular 
instabilities with ni — 3,4 (see Sect. 4.1). 



suits by Andersson (1998) and the results presented here 



4 ■4- Maximum and minimum mass of a neutron star 

Table |^ shows, that the maximum baryonic mass of a 
non-rotating (rotating) EPNS is 2.33M0 (2.45 - 2. SM©), 
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Table 5. Properties of the maximum gravitational mass configurations of non-rotating and with Kepler frequency 
rotating PNS's and NS's. 



EOS Mg Mb Rini Wc Tc Mg Mb Rint J 

[Mq] [Mq] [km] [fm-3] [MeV] [Mq] [Mq] [km] [fm-^] [s"^] [Mq km] 



2.33 12.59 1.09 40.4 2.18 2.45 16.93 1.08 7609 3.06 

2.33 11.83 1.10 40.7 2.22 2.50 15.44 1.09 8803 3.69 

2.33 11.91 1.03 81.6 2.30 2.54 15.90 0.97 8602 4.29 

2.37 11.29 1.09 92.0 2.33 2.62 15.43 0.95 9067 4.61 
2.33 10.82 1.14 41.8 2.28 2.56 15.38 0.93 9384 4.49 

2.38 10.70 1.18 45.6 2.31 2.65 14.43 0.93 9701 4.94 
2.27 10.29 1.12 2.28 2.54 14.31 0.88 9956 4.98 
2.41 11.24 1.12 120.6 2.34 2.68 15.47 0.92 9073 4.80 
2.38 10.48 1.20 62.1 2.28 2.66 14.14 0.92 9864 5.03 
2.41 9.83 1.25 2.30 2.70 13.74 0.95 10631 5.64 



EPNS?£;,4 


2.05 


EPNS^to4 


2.05 


LPNSfL04 


2.09 


LPNS^Los 


2.08 


LPNS^Yo4 


2.05 


LPNS^Yo3 


2.04 


LPNS^Olo4 


2.05 




2.08 




2.01 


CNS 


1 99 



whereas the maximum baryonic mass of a CNS is 2.41 Mq 
(2.7 Mq) (see also Fig. |^). Since the maximum bary- 
onic mass of a LPNS and HNS is larger than the maxi- 
mum baryonic mass of a EPNS, too, it is obvious that a 
PNS based on our nuclear EOS cannot collapse to a black 
hole during its Kclvin-Hclmholtz cooling phase if any fur- 
ther accretion is neglected (see Baumgarte et al. 1996|) . 
The situation changes if hyperons, meson condensation, 
or a quark-hadron phase transition are included (see, e.g.. 
Brown & Bethe 1994; Pandharipande et al. 1995; Prakash 



et al. 1995; Glendenning 1995). Another important point 
concerning the maximum mass of a CNS is due to the dif- 
ferent maximum masses of EPNS's and CNS's rotating at 
their Kepler frequencies. The EPNS can support a bary- 
onic mass of 2.45 — 2.5Mq, which is only slightly higher 
than the maximum baryonic mass of the non-rotating CNS 
(2.41Mq), but 0.2 — O.25M0 less than the maximum bary- 
onic mass of the CNS rotating at its Kepler frequency. 
Such supramassive CNS's enter unstable regions during 
their spin down evolution due to unstability against ax- 
isymmetric perturbations and may finally collapse to a 
bla ck hole. For a disc ussio n of these supramassi ve CNS's, 
see 



Cook et al. (19941) and ^algado et al. (19941) 



In first approximation, the gravitational mass of non- 
rotating PNS's and HNS's increases quadratically with the 
entropy per baryon, which is taken constant throughout 
the star (see Prakash et al. 1997| ): 

Mcis) - Mg{T = 0) (1 + + . . .) . (13) 

We obtain for A approximate values ~ 0.010 — 0.011 
(~ 0.005) for the cases without (with) trapped lepton 
number (see Table ^ . These values are in agreement with 
the values derived by prakash et al. (1997 ). In the inves- 
tigation of Goussard et al. (1997) the gravitational mass 



shows the opposite behaviour: The gravitational mass de- 
creases slightly with increasing entropy per baryon. This is 
probably caused by the smaller temperature dependence 
of the EOS of |Lattimer fc Swesty (1991] ), as it was pointed 
out in Sect. In the case of rotating stars the behaviour 
seems to be more complex since the value of A is nega- 
tive for HNS'i and LPNS^Yo4 and positive for HNS"^ and 
LPNSylo4- Another interesting point is, that the maxi- 
mum gravitational and baryonic mass of the non-rotating 
PNS's are determined by the EOS used in the core and 
do not depend on the EOS used in the envelope (compare 
the models EPNS|^£j,4, EPNS^^£j,4, and LPNS^Yo4 in Ta- 
ble I and Figs. ^ and |l|). In the case of rotating PNS's, 
the maximum gravitational mass of the EPNS is however 
smaller compared to the corresponding LPNS model (see 
Table || and Figs. ^ and This behaviour is caused by 
the higher Kepler frequency which can be supported by 
the LPNS models (see Sect. O). 



The minimum baryonic mass of our non-rotating 
EPNS sequences is in the range 0.95 — 1.24Mq. The mini- 
mum mass is increased in the case of Kepler rotating stars 
to 0.97 — 1.29A/q). If accretion is neglected the baryonic 
mass is conserved during the evolution of the EPNS to 
the CNS. We follow the evolution of the minimum mass 
EPNS|.^£i4 and EPNS^4£i4 model to a CNS via the LPNS 
and the HNS stage (see Table 1). The minimum mass of a 
CNS born in a supernova is therefore determined by the 
minimum mass model of the EPNS sequence. The above 
ranges of baryonic mass correspond to a lower limit of 
the gravitational CNS's mass in the range 0.89 — 1.13Mq. 
This is by a factor of ten larger than the minimum mass 
of the CNS sequence. This property of EPNS's was also 
recently found by Goussard et al. (1998). 
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EOS (Mb = O.95OM0). 



Table 6. Evolution of minimum mass star for the non-rotating EPNSf£^4 EOS {AIb = 1.238M0) and the EPNS"''''^ 



YL04 



EPNS|?£;,4 EPNS|.to4 



EOS Mg Rini Tc Mg Rini ric 

[Mq] [km] [fm^^] [MeV] [Mq] [km] [fm^^] [MoV] 



EPNS 


1.207 


45.61 


0.340 


17.3 


0.942 


44.49 


0.281 


15.1 


LPNSfL04 


1.199 


19.77 


0.308 


33.8 


0.926 


25.34 


0.219 


26.7 


LPNS^Los 


1.186 


17.49 


0.327 


37.1 


0.922 


20.38 


0.250 


30.7 


LPNS^Yo4 


1.185 


15.83 


0.384 


18.9 


0.926 


18.05 


0.307 


16.1 


LPNS^Yo3 


1.168 


14.24 


0.404 


20.3 


0.915 


15.31 


0.331 


17.6 




1.168 


16.28 


0.341 


43.2 


0.911 


18.26 


0.267 


36.5 




1.146 


13.18 


0.429 


23.7 


0.898 


13.75 


0.354 


20.5 


CNS 


1.134 


11.53 


0.481 





0.890 


11.58 


0.410 






Variations of the location of the transition region be- 
tween the envelope (high entropy per baryon) and the core 
(low entropy per baryon) of the EPNS's do not change 
the minimum mass considerably. Starting the transition 
region at lower densities will lead to smaller minimum 
masses of the EPNS's, but the lower the initial mass of 
the EPNS is, the higher is the entropy per baryon in the 
envelo pe a nd in the core (see e.g. Keil et al. 1996 ; Po: 
et al. l |999| ) . This effect drives the minimum mass back to 
higher values, so that our results of the minimum mass 



range is a good approximation (see also Goussard et al, 
1998|3 



Timmes et al. (1996) examine the most likely masses of 
NS's using the numerical data of Woosley & Weaver (1995) 



who simulated Type-II supernovae with progenitor stars 
in the mass range between 11 and 40 Mq. They obtain a 
lower limit of the NS mass, which depends on the mass 
and the composition of the progenitor, of 1.15 — 1.27 Mq. 
This lower limit is comparable to our results. 



4-5. Sensitivity of the results 

In Table |^, we compare the properties of our PNS and 
NS models for a fixed baryonic mass, M-q = 1.5 Mq. This 
canonical value corresponds to the measured gravitational 
masses, 1.35ib0.27 A/q, of neutron stars in binary systems 



(see Thorsett ct al. 1993; Van Kerkwijk et al. 1995) 



As can be inferred by comparing the LPNS models 

LPNSylo4(64-62) ' LPNSylo4(64-22) ' ^^'^ LP^^YL04(64-63) 



with the model LPNSylo4 (^^e Sect. 2.2), the location of 
the neutrino sphere has nearly no effect on the gravita- 
tional mass and the central density. However, the circum- 
ferential radius and the Kepler frequency vary by up to 
10% and 15%, respectively. 



The use of an isothermal, instead of an isentropic, 
EOS in the envelope of the HNS models HNS'^^^si 
HNS""""^^^ has only small effects on the properties of the 
HNS's (see Table ||). If thermal effects in the envelope are 
however neglected (models HNS'^^'^i and HNS™"^)^ the 
circumferential radius is reduced by ~ 10%. This yields 
to an increase of the Kepler frequency by ~ 20 %. Though 
the assumption of zero temperature in the envelope does 
not change the resulting mass and central densities, the er- 
ror made in the circumferential radius and in the Kepler 
frequency might be rather large. 

5. Discussion and conclusion 

The aim of this paper was the investigation of the proper- 
ties of rapidly rotating PNS's consisting of nuclear matter 
(n, p, e~ , /i~ ) under the influence of trapped neutrinos. We 
used a recently developed nuclear EOS for h ot, dense mat 



ter in the nuclear Thomas- Fermi approach ( Strobel et al 



1999). The nuclear EOS was extended in this paper to 



subnuclear densities and to different compositions of PNS 
and HNS matter (i.e. trapped neutrinos, constant entropy 
per baryon, . . . ) . Our results for subnuclear densities are 



1998 ) and Gondck et al. (1997 ). However, considerable de 



comparable to the results derived by Lattimer fc Swesty 
1991 ) use d in the investigatio ns of Goussard et al . (1997, 



viations occur at densities around and above nuclear mat- 



ter density, since the EOS of Lattimer fc Swesty (1991 ) 



shows a smaller temperature dependence. In general, this 
difference to our results has only a small impact on the 
properties of PNS's, but impacts on HNS's with canonical 
mass. The impact becomes considerable for more massive 
stars. 

An investigation of this kind should investigate the 
properties of PNS's and HNS's for the whole range of pos- 
sible masses. So far lepton concentration profiles and en- 
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tropy per baryon profiles were derived only for one fixed 



mass in simulations of PNS evolution (e.g. 


Burrows & Lat- 


timer " 


986 




Keil & Janka 1995 




Burrows et al. 1995). In 



our calculation, the extension to the whole range of masses 
was done by assuming that the lepton concentration and 
the entropy per baryon does not depend on the mass of 
the PNS. 

As a result wc find that the minimum gravitational 
mass of a NS is determined at the earliest stage of a PNS, 
so that the mass of a NS formed in a Type-II supernova is 
larger than 0.89 - 1.13 Mq, which confirms similar r esults 
of a recent investigation by Goussard et al. (1998 ). The 
exact lower limit of the NS mass depends on the used 
entropy per baryon in the EPNS model. The quoted mass 
range was obtained by using s = 4 and s = 5 as lower 
and upper limit of the entropy per baryon in the envelope 
of the EPNS ( Burrows et al. 1995 ). The minimum mass 
of EPNS's is approximately by a factor of ten larger than 
the minimum stable mass of CNS's. 

The maximum possible baryonic mass of a CNS is al- 
ways larger than the maximum possible baryonic mass of 
the PNS's, this means that once a PNS was formed, it 
cannot collapse into a black hole, if further accretion is 
neglected (see also Takatsuka 1995| ; Bombaci 1996 ). One 
exception is the case of the most massive stars rotating 
near or at its Kepler frequency, which possibly collapse 
to a black hole during their time evolution, for a discus- 

[1994|). 



sion, see 



: ICook ct al. (1994D and ^algado et al. (1994| ) . This 
statement holds for stars with a pure nucleonic/leptonic 
composition. If one includes hyperons and/or quarks, the 
maximum baryonic mass of a CNS decreases (e.g. Hubei 
et al. 1998; Balberg et al. 1999) and may be smaller than 
the maximum baryonic mass of a PNS. Then, a sufficiently 
massive PNS may collapse to a black hole during delop- 
tonization (Prakash et al. 1997). The maximum gravita- 
tional mass increases slightly with increasing entropy per 
baryon in contrast to the maximum baryonic mass (see 



also Prakash et al. 1997). As it was pointed out, Goussard 
et al. (1997) got a different result due to the use of the 
EOS derived by |Lattimer fc Swesty (1991| ). 

It turned out that the influence of rotation has sev- 
eral impacts on properties of PNS's and NS's. Whereas 
the minimum mass changes only slightly due to rotation, 
the effect on the maximum mass is rather large, particu- 
larly for CNS's. The central baryon density is nearly un- 
affected by rotation at the early stages of the evolution, 
whereas the impact on the the later stages is rather large. 
The effect on the circumferential radius decreases with in- 
creasing mass and the impact on EPNS's and LPNS's is 
slightly larger than the impact on CNS's. 

Wc investigated the influence of different shapes of the 
neutrino sphere on the structure of LPNS's. As expected, 
we obtained considerable differences (up to 10%) in the 
circumferential radii and, due to this, in the Kepler fre- 
quency. Other properties, as the central density, are barely 
changed by the location of the neutrino sphere. We have 



also considered different temperatures in the envelope of 
LPNS's and HNS's. If a very high temperature of 0.6 MeV 
is used in the envelope instead of a constant entropy per 
baryon, the gross properties of the LPNS change by less 
than 1 %. Larger deviations were obtained if the enve- 
lope is assumed to be cold. The flnite temperature effects 
should therefore not be neglected in the envelope. Further- 
more, we obtain that the thermal effects are comparable 
to the effects due to trapped neutrinos. This is in con- 
trast to the results of Prakash et al. (1997), who found 



the thermal effects to be much smaller than the effect of 
high lepton numbers. 

In our model we assumed that accretion onto the 
emerging NS stopped after the formation of the EPNS 



(see Burrows et al. 1995). Furthermore, we kept the angu- 
lar momentum constant during the deleptonization and 
the thermal cooling period. Under these presumptions, 
we obtain a lower limit on the periods of young NS's 
with a typical baryonic mass of 1.5 Mq between 1.56 
and 2.22 ms, which is in accordance with similar values 
obtained by Goussard et al. (1997 , 199^ ). These results 
support strongly the hypothesis that millisecond pulsars 
were accelerated due to accretion. With the same reason- 
ing, one obtains an upper limit of the stability parameter 
/3 < 0.082, which is smaller, for almost all NS's, than 
the critical value for the onset of dynamical and secular 
instabilities. Kepler rotating massive CNS's possess sta- 
bility parameter values up to /3 ~ 0.13 and might there- 
fore be secular unstable against m — 2 and to = 3 non- 
axisymmetric perturbations. 

Though we have assumed uniform rotation, it seems 
to be very reasonable that PNS's and young NS's rotate 
differentially (e.g. Janka fc Monchmeyer 1989|). As i t was 
found by [Schaab (1998D and [Goussard et al. (1998| ), dif- 
ferential rotation may considerably effect the structure of 
PNS's and NS's, and we will pursue our investigations 
to differential rotation in the future. Another issue that 
will be addressed in a future work is the effect of addi- 
tional degrees of freedom (e.g. hyperons) on the evolution 
of PNS's. 
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Appendix A: Nuclear equation of state 



The momentum and density dependent interaction is given by [the upper (lower) sign corresponds to nucleons with 
equal (unequal) isospin] ( [Myers fc Swiatecki 199C ): 



Vl2r 



no 



1 



1 



5(n2)U(lT0«-:^(lTC) /3 — 



Pl2 



PO 



Po 



\P1 



(A.l) 



The quantities no, po, and Tqt{~ Po/'^n^r) [To — p^/ifa with ffi — 0.5(rn,i -t-mp), see (A. 4)] denote the baryon number 
density, Fermi momentum, and the kinetic single-particle energy of symmetric nuclear matter at saturation (r denotes 
the isospin), respectively. The choice ^ 7^ C leads to a better description of asymmetric nuclear systems, and the 
behaviour of the optical potentials is improved by the term (T(2n/no)^/'^, where n^/^ = 0.5(nJ^'^ -I- ri^^^)- 

The parameter set is given by the most recent adjustment of Myers fc Swiatecki (1996 , 199^ ): a = 1.94684, 
(3 = 0.15311, 7 — 1.13672, a = 1.05, ^ — 0.27976, and C — 0.55665, which leads to the following properties of symmetric 
nuclear matter at saturation density (no = 0.16114 fm~^): energy per baryon uq = —16.24 MeV, incompressibility 
K = 234 MeV, symmetry energy J = 32.65 MeV and effective nucleon mass m*/mr — 0.867. 

The potentials radial dependence, g, is chosen to be of Yukawa type. The function is normalized to unity: 



+00 



i(ri2)dVi = 1. 



This leads to the following single particle energy u ( strobel et al. 1999 ) 

2 



n(27r)3 

with the one-particle potential Vt{pi). 

; 2 2To, 

{2Try no 

2 2To '■+°° 
' (27r)3^ 



E 



9 

Pi 

2m 



; + ^^r(pi))/r(pi)dV, 



(A.2) 



(A.3) 



Oil - A 



P12 

Po 



au - /3u f — ) 



■71 



7u 



Po 

\P12\ 

Po 

\P12\ 



2n 

no 

2n\ 3 



/r(P2)d^P2 



— /_,(p2)d3p2, 

no/ 



(A.4) 



[t and — T denote different isospin and the following abbreviations were used: a\ = 0.5(1 — ^)q;, f3\ = 0.5(1 — 0/3, 
71 = 0.5(1 -C)7, (Ji = 0.5(1 -C)cr, au = 0.5(1 +0", Pu = 0.5(1 + C)P, 7u = 0.5(1 + C)7 and a,, = 0.5(1 + C)cr] and the 
baryon number density: 



(2^)^ 



/+00 
/r(pi)d^Pl. 



(A.5) 



In the Eqs. (A.3), (A.4), and (A.5) /t(p) denotes the Fermi-Dirac distribution function (see Fig. A.l as an example) 
of a baryon with isospin t (fee = 1): 



Mpi) 

where e,- denotes the one-particle energy: 



1 + exp 



(l(e.(pi)-M;)))"\ 



er(pi) 



A. 

2mr 



VApi). 



(A.6) 



(A.7) 



Hint: /i^ in Eq. ( A6) is not th e chemical potential in no rmal sense, because of the density dependent part in the 
interaction, for an explanation see Myers fc Swiatecki (1990 ), Appendix A. The chemical potentials of neutrons and 
protons, /ir, can be derived over the thermodynamic derivatives: 



) +l]u + mr = \n[ ) +1] f + nir 



d 
drir 



drir 



(A. 



in which / is the free energy per baryon and to^ the rest-mass of neutrons or protons. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

en/e„(T=0) 

Fig. A.l. Fermi-Dirac distribution function for the neu- 
trons of the LPNS^\o4 EOS at n = 1 fm--\ T = 79.8 
MeV, Tf = 278.2 MeV, T/Tp « 0.29, n'^ = 512.8 MeV 
and en(T = 0) = 521.2 MeV at the Fermi surface. 
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